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Abstract 



In principle, observables as for example the sphaleron rate or the 
tunneling rate in a first-order phase transition are gauge-independent. 
However, in practice a gauge dependence is introduced in explicit per- 
turbative calculations due to the breakdown of the gradient expansion 
of the effective action in the symmetric phase. We exemplify the sit- 
uation using the effective potential of the Abelian Higgs model in the 
general renormalizable gauge. Still, we find that the quantitative de- 
pendence on the gauge choice is small for gauges that are consistent 
with the perturbative expansion. 
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1 Introduction 



Functional methods are an indispensable tool in studying the physics of spon- 
taneous symmetry breaking and phase transitions [TJ [2j [3J, H] - The main 
advantage in this approach is that the effective action encodes the ground 
state of the system in a transparent way. On the technical side, it facilitates 
the resummation of tadpole diagrams and reduces perturbative calculation 
to the subset of one-particle-irreducible diagrams. 

A large drawback of the effective action is however that it is not explicitly 
gauge- independent El [TJ El EJ [10]. This makes it necessary to distinguish 
between the gauge dependence that is expected from the one that is intro- 
duced by approximation schemes that break the gauge invariance addition- 
ally. The first class of gauge dependences is well represented by the Nielsen 
identities while the second can for example arise from the use of a loop ex- 
pansion in perturbation theory or from the expansion of the effective action 
in gradients. 

The aim of the present work is to disentangle these different sources of 
gauge dependence for vacuum transitions [TTJ at finite temperature in the 
Abelian Higgs model. After a general introduction to the model (Section 
[2]) and the Nielsen identities (Sections [3] and @J, the effective potential is 
calculated (Section [S]). The main focus is hereby on resummation of infrared 
effects at finite temperature. In this section, it is also explicitly demonstrated 
that the position of the minimum of the effective potential transforms under 
changes of the gauge fixing parameter £ according to the Nielsen identity. 
Subsequently (Section^, the same is demonstrated for the (off-shell) effective 
action in the gradient expansion. Finally (Sections [7] and E]), the gauge 
dependence of the tunneling action and the sphaleron energy are discussed. 
In these cases an additional complication arises, namely the breakdown of 
the gradient expansion in the symmetric phase. We discuss, to what extend 
the gauge dependence of these phenomena can be quantified and potentially 
reduced before we conclude (Section [9]) . 

2 The model 

In the following we lay out the details of the model under consideration. 
In order to simplify the analysis and to focus on the main impact of the 
gauge dependence of the effective action, we discuss an Abelian Higgs model. 



2 



This model has all necessary ingredients that occur in the Standard Model 
but does not contain fermions or the correct symmetry breaking pattern as 
observed in Nature. Nevertheless, for the order of perturbation theory we 
work at these features are not important and the presented arguments can 
immediately be carried over to the Standard Model case. 
The Lagrangian is given by 

C = D^*D^ - -F^F^ - , (1) 

with the scalar potential 

^($*$) = -^2$*$ + A($*$)2 = ^($*$ _ v 2f. ( 2 ) 

The gauge covariant derivative is given by = (<9 M — igA^)^, while the 
field strength tensor reads F^ v = d^A v — d v A^. 

We are interested in the effective action Y that is obtained [H |3| by 
Legendre transformation of the generating functional of connected Greens 
functions W(j) 



(3) 

<f>=dW/dj 



e iW(j) _ J D^DAe iSdixCm+jq> 

r(0) = w{j)- j d 4 xj^ 

This implies after a shift in the integration variable $—?■$ + <fi 

e ir = J V$VA e { S di *^+*)+i * . (4) 

By construction the field in the shifted theory does not have a vacuum expec- 
tation value (vev). Hence, the bi-linear term containing the source j has to 
cancels all tadpole diagrams and beyond this cancellation, the source j has 
no impact on the effective action. Therefore V contains only the connected 
one-particle-irreducible vacuum diagrams in the shifted theory. The effective 
potential is obtained by restricting the effective action to homogeneous field 
expectation values (ft, V e ff = Y / J d A x. 

For perturbative calculations, the Lagrange density has to be supple- 
mented by a gauge fixing procedure that results in an additional gauge fixing 
term and a contribution from the Fadeev-Popov ghosts. 
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A common choice is to use the general Lorentz gauge 

Lorentz gauge : £ gf = - — (d^A^) 2 , 

Cpp = —cDc. (5) 

In this gauge the Fadeev-Popov ghosts do not couple to the remaining par- 
ticles, but the Goldstone boson and the unphysical k^ polarization of the 
gauge boson mix after spontaneous symmetry breaking. Nevertheless, this 
problem is removed in the Landau gauge, £ — > 0. 
Another common choice is the i?g-gauge [12] 

i^-gauge: C gf = ~{d^ + z<?£(<f$ - $ *0)) 2 , 

C FP = -c(D + ^ 2 (0*$ + $*0))c. (6) 

The contribution containing the background field <fi is designed to cancel 
the mixing term between the Goldstone mode and the polarization of 
the gauge boson thus partially diagonalizing the propagator structure of the 
theory. However, having an explicit dependence on the background field in 
the gauge fixing term leads to some complications. 

First of all, in Lorentz-gauge a shift in <fi is equivalent to attaching an ex- 
ternal $ field to the diagram [T3l H] (this follows immediately from inspecting 
the Feynman rules). Hence the tadpole can be obtained from the derivative 
of the effective action. In general gauges this only holds in the minimum of 
the potential (where the tadpole has to vanish) while in Lorentz gauge this 
is also true away from the extrema of the potential. 

One consequence of this fact is that the two-point functions of the Higgs 
and Goldstone bosons (at vanishing external momentum) can be obtained 
from the second derivatives of the effective potential. In particular, the 
Goldstone bosons are massless in the broken phase (what can be a curse 
or a blessing). In _R^-gauge the Goldstone bosons have a mass already on 
tree-level for £ ^ 0. 

Furthermore, any tadpole that is generated in a higher loop perturbative 
calculation can in Lorentz gauge be absorbed by adjusting the appropriate 
counterterms. This is particularly handy in zero temperature calculations 
where the Higgs vev is an input observable parameter and kept fixed order 
by order in the perturbative expansion. Especially when several background 
fields are present, it is not guaranteed in the .R^-gauge that the tadpoles 
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can be absorbed into the counterterms of the Lagrangian what makes it 
technically more difficult to keep the vevs fixed in higher loop calculations 
[13]. Note that at finite temperature this feature is not so important, since 
part of the shift of the vev by loop effects is physical and one has to deal 
with a temperature-dependent vev anyway. 

The behavior of the effective action under gauge transformations was 
originally analyzed in the Lorentz gauge by Nielsen [Sj. Similar relations 
have been derived shortly after for i?^-gauges P [7J IE] but their application 
is somewhat more involved. Still the general picture is the same as in the 
case of the original Nielsen identities: While the value of the effective action 
in its extrema is gauge-independent, the position of the extrema can have a 
gauge dependence. This implies that the critical temperature of the phase 
transition is gauge-independent. For vacuum transitions, the configurations 
that extremize the effective action can be gauge-dependent while the value 
of the action for the configuration is not. The prove of this is sketched in the 
next section. 

Ultimately, we aim to discuss the gauge dependence of the sphaleron 
energy and the tunneling rate. In this work we focus on the i?^-gauge for two 
reasons. First, we find it technically easier to eliminate the mixing between 
the gauge boson and the Goldstone boson and to deal with the peculiarities of 
the i^-gauge than to perform loop calculations in the general Lorentz gauge 
0. Second, the convergence of the perturbation theory is in the Lorentz gauge 
less obvious than in the .R^-gauge [TU [TJ)J [T6] . 



3 Nielsen identity in gauge 

In this section we briefly review the Nielsen identities for the effective poten- 
tial, paying special attention to the additional complications arising in the 
i?g-gauge. In order to derive the Nielsen identity for a general gauge fixing 
we write the latter in the form, 



2i 

c FP = -c(ga M + ^ + ^ H ^)] c , (7) 



1 See [TT| for a zero temperature analysis of the gauge dependence of the effective action 
in general Lorentz gauge. 
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where F = <9 M A M in Lorentz gauge and F = d^A^ + ig£((f)*$ — $*0) for 
the .R^-gauge. For the derivation of the Nielsen identities, it is convenient 
to discriminate the field expectation value = (<&) that appears due to 
the Legendre transformation of the generating functional from the explicit 
dependence on the background field in the gauge fixing term, which is denoted 
by here and can be considered as an external parameter at this stage of 
the calculation. Both fields will be identified in the end, but for now we 
distinguish them. 

The dependence on F in general leads to a gauge dependence of the effec- 
tive action V . For an arbitrary change in the gauge fixing term parametrized 
by F — > F + SF and £ — > £ + the change of the functional W is given 
by 

SW= -iJd A xJ d* yjl (8 g <S> t (x)c(x)c(y)8'F(y)}, (8) 

where 5'F = SF — F/(2£)8£. Furthermore, we collectively denote by 0j all 
field expectation values that the effective action depends on, and by 5 g §i the 
variation of the corresponding field operator under a gauge transformation 

5 g $ = ig$, 5 g ^* = -ig$; 8 g A^ = «9 M . (9) 

The expectation value of an operator (O) is defined in the usual path integral 
sense. 

The functional W is gauge invariant up to the source term involving ji 
and the gauge fixing term C g f + C FP . The identity (JS} reflects the fact that a 
change in the gauge fixing functional can via a gauge transformation shifted 
into a change of the source term. Now, notice that for any external parameter 
A one has the relation 



dX 



dW{j,X) 



<f>=const ^ 



(10) 

j=const 



what can be verified using the definition ([3]). The gauge fixing parameters £ 
and are external parameters such that one can translate (jHJ) into a relation 
for T if the left hand side is understood as being varied with kept fixed, 
namely 

5T = i Jd 4 xJ d A yj^(5 g ^Hx)c(y)5'F(y)) . (11) 
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This immediately yields the Nielsen identity for the effective action, which 
results when only considering a change in the gauge parameter £, 



«f = " / A {-^) CMx) + 1W) C * (X) + W^) CAx) } • (12) 

where the coefficients, taking a possible £ dependence of F into account, are 
given by 

C A ,(x) = l -jd'y(d,c(x)c(y)(F(y)-2^dF(y)/dO), 



i 



C+(x) = ~ I d 4 y{ig®(x)c(x)c(y)(F(y) - 2£dF(y)/d£,)) 



C r (x) = % -J d k y((-ig^(x))c(x)S(y)(F(y)-2^dF(y)/dC)). (13) 

For a constant field expectation value, chosen to lie along the real axis <fi = <p*, 
and vanishing background gauge field = 0, one obtains the well-known 
Nielsen identity for the effective potential, 

5 M + C «!.o. (14) 

The coefficient C is obtained when evaluating C{x) = + C^*\^=^* for a 
constant field expectation value <fi. In i?^-gauge, the general expression for 
C{x) reads 

C(x) = J d ±y (c(x) X (x)c(y) (dpAHy) + y/2g^(v)) ) , (15) 

while in Lorentz gauge the second summand is absent. Here we have inserted 
the decomposition $ = <p + {h + z'x)/V2 where h and x denote the Higgs 
and Goldstone field operators, respectively. Note that the restriction to one 
real background field and vanishing gauge fields can lead to spurious minima 
in the effective potential but is justified in our class of gauges [5]. This is 
mostly due to the invariance of the full Lagrangian (including gauge fixing 
and sources) under the transformation $ — > $* and — > —A^. 

For the Lorentz gauge these identities can be used directly when inserting 
the gauge fixing term F = d^A^ 1 . In the -R^-gauge, the effective action 
depends in addition parametrically on the external field 0, i.e. V = r[0;</>]. 



7 



We are ultimately interested in deriving a Nielsen identity for the effective 
action F[0] = r[0;0]|^ = .. Compared to the Lorentz gauge the complication 
arises that the field derivatives appearing in Eq. ([12]) act only on <p, while the 
field derivative of F[0] contains also a part proportional to the derivative of 
the effective action with respect to the background field <fi. Nevertheless, it is 
possible to derive a Nielsen identity also for F [</>]. The key observation is that 
the dependence on the background field enters only via the gauge fixing 
term [8]. For that reason, the change of the effective action when varying the 
background field also obeys a Nielsen identity, which can be directly obtained 
from Eq. flTT), 

m =i f^)( SMx * xmB -w)- (16) 

An analogous relation holds for the derivative with respect to the complex 
conjugated field. One can use the relation above to express the derivatives of 
F in terms of the derivatives of T with respect to the field expectation values 
0, 0* and A M , 

6T f .„ x-^ „ , . ST 



my) S ^ ^J isM W^^' (17) 



where 



(^(y,x) = S ij S(x-y) + i(s g 9 j {x)c(x)c{y)?^^ . (18) 

These relations imply that, if the fields fulfill the equation of motions de- 
rived from the effective action r[0;0], i.e. 5T/5(pj = 0, then the same field 
configurations also correspond to a stationary point of the effective action 
F [</>]. In addition, if C is invertible, it follows that the effective action F[0] 
also fulfills a Nielsen identity which is of the same form as Eq. f[T2"j) except 
that the coefficients are replaced by the coefficients (7j given by 

Ci(y)= [ d*xC 3 (x)(C-%(x,y), (19) 



for i = A^, <f>, (ft* . Concretely, for the case of the i?^-gauge considered here, 
and assuming a vanishing background gauge field, one has 



dj = Si^x-y) (20) 

mx)c(x)c(y)<$>*(y)) -($*(x)c{x)c(y)$*(y)) \ 
-($(x)c(x)c(y)$(y)) (<S>*(x)c(x)c(y)<S>(y)) . J 



+ ^g 2 ^ 
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Figure 1: One- loop diagram contributing to the Nielsen coefficient Co- 

When considering an expectation value that is constant in space-time, the 
right-hand side can depend only on the difference of the coordinates, x — y, 
and the equations can be simplified. The resulting Nielsen identity for the 
effective potential for constant and real field expectation value in i?^-gauge 
has then the form 

^ + Q,M = „, (21) 

where V eff (<f>) = K//0, 4>)\^ =<t> , and 

C, = cJl-ig 2 i [ d\x-y)(x(x)c(x)c(y) X (y))] ■ (22) 



The leading contribution arises at one-loop level. For a constant background 
field, only a single diagram contributes to C in the .Rg-gauge, which can be 
represented by the diagram shown in Fig[U The coefficients Co and Co differ 
only at higher orders, which we do not consider in the following. Therefore, 
we use the notation C also in the case of .Rg-gauge. Before checking ( |2T!) 
explicitly in Section [5j we give attention to an analogous identity for the 
kinetic term in the next section. 



4 The Gradient Expansion of the Nielsen Iden- 
tity 

As we have just discussed, the effective action contains a certain gauge de- 
pendence. Nevertheless, the action evaluated in a saddle point is gauge- 
independent what is for homogeneous fields reflected by the Nielsen identi- 
ties. This is explicitly demonstrated in the subsequent section for the broken 
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phase of the effective potential where the convergence of perturbation the- 
ory is well under control. However, the sphaleron energy is not a spatially 
constant configuration and in particular also depends on the effective action 
evaluated in the symmetric phase. This is problematic, since even though 
the sphaleron energy is gauge-independent in principle, a gauge dependence 
can be introduced due to the breakdown of the gradient expansion. Still, 
even for configurations that do not fulfill the equations of motion, a gauge- 
transformation acts on the (off-shell) effective action according to the Nielsen 
identity. In the following we derive this relation for the effective action in 
gradient expansion following [TT] . 

To be specific, consider the effective action in gradient expansion 

d A x (Z((f>) d^4> - V eff {<f>) + 0(d 4 )) , (23) 

where we introduced a renormalization Z of the kinetic term that in the 
sphaleron case depends on temperature T, the field value and the gauge 
parameter £. The coefficient C(x) from the last section that enters the Nielsen 
identity for the effective action depends functionally on the field value. It 
can be expanded in field gradients as [TT] 

C(x) = C o (0) + D{4>)d^4> - d^D(<p)d^] + 0(d A ) , (24) 

where the coefficients depend on the field value evaluated at position x. Note 
that the term involving D corresponds to a total derivative. Nevertheless it 
gives rise to a non-zero contribution to the right-hand side of the Nielsen 
identity for the effective action, and it turns out that it has to be taken into 
account in a consistent gradient expansion. This term has been neglected 
in [11] but was not important there at leading order due to the different 
counting A ~ g 4 . The Nielsen identity is fulfilled when in leading order in 
gradients the identity (121 j) for the effective potential is established, while in 
the next order in gradients one finda^] 

Notice that this relation cannot ensure the gauge independence of the 
tunnel action or the sphaleron rate, since the gradient expansion does not 



2 At finite temperature the correction to the kinetic term can differ in general for the 
temporal and spatial components. The Nielsen identity can be easily generalized to this 
case. In the following, Z refers to the correction corresponding to the spatial components. 
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apply in these cases, due to d^d^ifi ~ V e ff- Still, this relation constitutes an 
essential check of the consistency of the perturbative scheme that is used to 
evaluate the effective action. The topic of the subsequent sections is to study 
this identity in the Abelian Higgs model close to the broken phase explicitly, 
where perturbative evaluation of all quantities is plausible. 



5 The effective potential to order g 3 and A 

In the following we reproduce the effective potential to order g 3 and A close 
to the broken phase. Subsequently, we discuss its gauge dependence and the 
corresponding Nielsen identity. Before we do so, we briefly motivate why a 
counting g 3 ~ A is the appropriate choice in the context of cosmology and 
finite temperature. It is well known that in the present model the strength of 
the phase transition crucially depends on cubic contributions to the effective 
potential of the form (g<f)) 3 T. At the same time the strength of the phase 
transition is for viable baryogenesis [TO Q2] constrained by 4> C /T ~ g 3 /\ > 
1. So the largest value of A that is interesting in cosmology is of order 
g 3 . Besides, for larger values of A, the convergence of perturbation theory 
becomes worse and completely breaks down at A ~ g 2 . For smaller values 
of A, the convergence of the perturbative expansion improves. However, the 
high temperature expansion that we want to employ in the following would 
break down since for A ~ g A , we find my/ /T ~ g<p c /T ~ 1. Therefore, we 
count in the following <p c ~ T and A ~ g 3 as was done in pj9j %. 
At one-loop order the effective potential is given by 

VeffW = V{4>) + V ni lB(m 2 /T 2 ) , (26) 

i 

where the sum runs over all species, ni depends on the statistic of the field, 
and the function Ib is given by 

Ib(v) = I dx %2 log ( 47f2n2 + x2 + v) 

n J 

- const + ±-y - ^y 3 ' 2 + 0{y 2 ) . (27) 

12 67T 

Notice that the zero temperature contributions are of order y 2 and are ne- 
glected in the following. 



3 At zero temperature, the counting A ~ <? 4 is more appropriate fTj. 
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The tree level inverse propagators of the Higgs, Goldstone, gauge fields 
and ghosts are respectively 

^h/x/FP = P 2 ~ m h/ x /FP J 

iP A X = p 2 g^-( 1 - 1 /0p»Pv- m A9^, 

where we have introduced the field-dependent masses 

ml = ^ 2 -v 2 ), 

m\ = ^(0 2 -^ 2 )-2^ 2 2 , 

m\ = 2g 2 <p 2 , 
m 2 FP = 2^ 2 2 . (28) 

These expressions can be used to write the effective potential in the mean 
field approximation (using the leading terms in ( )27|) ) 

^-v 2 ) + 2tf<p 2 

+ 2[3/0 2 + ^0 2 ] 
= ^(0 2 -. 2 ) 2 + ^(A(20 2 -, 2 ) + 6,V), (29) 

where the different contributions are from the Higgs, the Goldstone, the 
gauge bosons and the ghosts, respectively. 

Notice that at tree level and for = v , the poles of the Goldstone, 
the polarization and the ghosts are all situated at p 2 = 2£g 2 2 and the 
corresponding contributions to the effective action cancel @. However, at finite 
temperature the minimum of the effective action moves away from <p = v and 
higher loop corrections become important that have to be resummed. In the 
following we show that these two effects indeed cancel each other in leading 
order. 

The necessity of resummation [21 121] at finite temperature arises because 
of diagrams as depicted in Fig. |5J In general, all finite temperature contri- 

4 A term that depends on £ but not on <j> can be canceled by an appropriate measure 
for the ghost fields [8]. 



V mean - f , eld = \(tf-v 2 ) 2 + ^^ 2 -v 2 ) + 
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Figure 2: The daisy diagrams that are resummed. 
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Figure 3: The leading contributions to the self-energy of the gauge bosons. 

butions to the self-energies are UV finite. Once the sum over the Matsubara 
frequencies is performed (or if the real time formalism is used), the integrand 
contains the particle distribution functions that are exponentially suppressed 
for momenta larger than the temperature. Hence, the graphs that are ap- 
parently UV divergent can be estimated to be of order of the temperature. 
In particular, tadpole diagrams of the self-energy that arise from the gauge 
interaction are of order g 2 T 2 (e.g. the contributions to the self-energy of the 
gauge bosons shown in Fig. [3]). 

If the particle in the loop has a mass mi and the self-energy is of order 
g 2 T 2 , adding self-energies leads to additional factors 



As long as n > 0, this yields only a subleading correction of order g 2 . Still, 




(30) 



(27T77.T) 2 + p 2 + m 2 
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for the zero mode n = 0, this contribution can be sizable and the according 
diagrams have to be resummed. This problem is particularly severe when 
mf <C g 2 T 2 and the convergence of the perturbative expansion crucially de- 
pends on this resummation. This is for example the case for the longitudinal 
gauge bosons, that in the symmetric phase are massless on tree level, but 
receive self-energy corrections of order S ~ g 2 T 2 at one-loop. But even in 
the broken phase the self-energy and the mass are both of order g 2 (p 2 and 
resummation is essential. 

In order to resum these diagrams consistently, one has to absorb at least 
a part So of the self-energy into the propagator. In the broken phase, the 
details how this is done do not really matter as long as 

ml ff = m} + E >E-E . (31) 

This ensures that IR divergences are tamed and the perturbative series con- 
verges. These self-energy contributions depend on the temperature and they 
have to be compensated by introducing a counterterm of equal size in order 
to avoid a temperature-dependent regularization scheme. These countert- 
erms contribute at higher loop level as we will see below. So can in principle 
depend on momentum and also on the Matsubara number n. A particu- 
larly simple choice is to only absorb the self-energy of the zero mode in the 
propagator [TU] . 

For the choice of resummation scheme, it is crucial that in the present 
context we are only interested in static quantities. If non-static quantities 
as the plasmon damping rate are considered, the resummation of the zero 
modes is not enough to ensure the convergence of perturbation theory as 
discussed in [22] • In this case more elaborate schemes like hard thermal loop 
resummation have to be employed [23J. 

In the broken phase (and for £ ~ 1), it suffices to absorb the leading term 
of the self-energy of the gauge bosons. After this is done, the pole structure 
of all gauge boson propagators is proper in leading order and corrections 
can be treated perturbatively. On the other hand, the symmetric phase is 
more problematic because the transverse polarizations remain massless in 
perturbative calculations. 

Within the context of non-Abelian gauge theories, adding additional 
lines using the self-interaction vertex of the gauge bosons results in a fac- 
tor g 2 T /rriA what is of order g in the broken phase. In the symmetric phase, 
on the other hand, the mass of the transverse polarizations is known [23| [24] 
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to be non-perturbative and only of order g 2 T what leads to the well known 
Linde's problem. For the Abelian case, there is no self-interaction of gauge 
bosons and the mass of the transverse polarization should vanish in the sym- 
metric phase (25]. But also here the convergence of the perturbative expan- 
sion is not obvious [261 121] • 

Compared to the unresummed one-loop result in Eq. ( 12"7|) resumming only 
the zero modes gives an additional contribution 

f d 3 p p 2 + m 2 + S 

T / 7T~T3 lo S 2~i 2 — • 32 

J (Ztt) 6 p l + m l 

One way to evaluate this expression is to determine its derivative with respect 
to m 2 which leads to 



T 



d 3 p ( 1 1 



(27r) 3 \p 2 + m 2 + S p 2 + m 2 
= -— (V™ 2 + S - V^ 2 ) ■ (33) 

Therefore, the resummation of the zero mode only affects the cubic term in 
Eq. (127]) . The resummed gauge boson masses are given by 

m% = 2g 2 <f) 2 , 

m 2 L = 2g 2 4> 2 + a 2 g 2 ct> 2 . (34) 

Here, we parametrized the thermal mass of the longitudinal gauge boson with 
the parameter a because it depends e.g. on the fermionic matter content of 
the model and is quite different in the Standard Model than in our example 
calculation. 

Now consider the Goldstone bosons that seem to be more interesting than 
the gauge fields due to the £ dependence in their mass. The Goldstone bosons 
are massless on tree level for £ = 0, such that resummation can have a large 
impact on them. If the leading high temperature contributions are taken into 
account, the self-energy arises from the diagrams in Fig. 2 and the effective 
mass is given by 

orders 2 : m 2 ;6// = ^(0 2 - v 2 ) + 2frY + ^(A + 3<? 2 ). (35) 



This fits nicely with the mean field result in Eq. ([29]) . If our counting was 
A ~ g 2 , the mean field result would give the correct leading order result for 
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Figure 4: The leading contributions to the self-energy of the Goldstone bosons. 

the effective potential and the above expression shows that the Goldstone 
boson has a mass 2£g 2 2 to this order. This ensures that the cubic terms 
of the Goldstone boson and the ghost cancel in the broken phase once the 
next to leading order is taken into account. However, our counting is A ~ g 3 , 
and to leading order the cubic terms coming from the three physical gauge 
bosons are important. Even though this term is ^-independent, it induces 
an additional shift in the vev such that the masses of the Goldstone boson 
and the ghosts do not coincide any more in the broken minimum of the 
potential. This problem is particularly severe for small £, since the ghost 
becomes massless while the Goldstone boson has a mass of order 

so the cubic contribution from the Goldstone boson is (partially) screened 
while the contribution from the remaining ghost is not. This introduces a 
gauge dependence of order g 3 into the effective potential. 

The solution to this dilemma is that there is a subleading contribution 
from the gauge bosons to the self-energy of the Goldstone bosons. This is 
easily seen by realizing that the integral that arises in the tadpole contribu- 
tion 

1 . f 1 

My) = tt^ Y] / dxx2 A 2 2 i — 2~— ' ( 37 ) 
^ ^ n + x + y 

is nothing else than the derivative of the one-loop vacuum contribution in 
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Figure 5: Two- loop contribution to the effective potential. The box denotes the 
counterterm of the thermal mass. 

Eq. (J2ZD 

Mv) = jjjlaiv) ■ (38) 



Hence there is a contribution from the two-loop diagrams [19J depicted in 
Fig. [5] that is of order g 2 rriAT 3 ~ g 3 T 3 cf). If £ < g this contribution needs 
to be resummed. This then ensures that the Goldstone boson has the same 
mass as the ghost in the minimum of the potential to order g 3 and their 
cubic contributions cancel each other once they are taken into account. The 
resummed Goldstone mass is given by 



\ T 2 
Tg 2 



orders 3 : m^ e// = -(0 2 - v 2 ) + 2^ + ^(A + 3g 2 ) (39) 



An 



(2m T + m L ) . 



In the following, we denote this expression for the Goldstone mass by m 2 . 

Finally, we comment on the contribution from the Higgs. Because the 
Higgs has a mass of order A0 2 , its cubic contribution and thermal mass are 
not relevant to our analysis. Still, there is a subtlety coming from the fact 
that the loop contribution to the self-energy in Fig. [6] is of order \ 2 <ft 2 T/m x 
and diverges in the limit £ — > 0. This divergence is due to the fact that the 
self-energy is evaluated at vanishing external momentum. If it is evaluated 
with external momenta close to the Higgs mass, it is of order X 2 (p 2 T/k ~ 
\ 3 l 2 T(f) what is small compared to the tree level mass and no resummation is 
necessary. Anyway, the numerical impact of this contribution is very small 
as long as £ is not strictly zero. 
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Figure 6: The Goldstone contributions to the self-energy of the Higgs boson. 

In conclusion, we arrive to the following expression for the effective po- 
tential to order g 3 and A 

Vto»,A) = ^(0 2 -, 2 ) 2 + g(A(20 2 -, 2 ) + 6,V) 
T 

(2m% + m\ + m 3 - m 3 FP ) . (40) 



12tt 



We emphasize that the resummed mass of the Goldstone boson is given by 
(|39|) . As argued before, the contribution from the Goldstone boson plus ghost 
is small when evaluated close to the minimum of the potential. For small 
£ <C X/g 2 it is of order A 3 / 2 , while for £ ~ 1 one finds 

sv ~ (K) 3 / 2 - (m| P ) 3 / 2 ) ~ -^f K - ^p) > ( 41 ) 

which is of order g\ ~ g 4 and also subleading. However, if it is not neglected, 
it induces a slight shift in the minimum of order 

'^^I^- (42) 

A shift in the position of the minimum is expected on general grounds, be- 
cause the field value is not a physical observable. Note that close to the 
critical temperature this shift is sub-leading, 



9 



V 7 ?, (43) 



because in our counting 4> C /T C ~ 1. Other contributions to the potential of 
order gX could eventually reduce this shift but actually this shift persists 
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even if higher order contributions are taken into account. This can be seen 
by inspecting all two-loop diagrams. It turns out that even though some 
contributions are nominally of order g A , they either do not depend on £ 
or are not linear in the temperature, so formally they cannot cancel ( l4Tj) 
completely. However, one expects at higher order a contribution of the form 

^"-It^< < 44 > 

to the potential close to the broken phase that removes the remaining gauge 
dependence of the effective action evaluated in the minimum. That this 
term is not included in the potential (HOI) introduces a subleading gauge 
dependence in the critical temperature of order 

T ~ 256vr 2 ' 1 j 

what is nominally of order g 3 and numerically very small. 

Let us now discuss to what extent the potential (|40j) respects the Nielsen 
identity (JHJ). First, we show that the shift (142]) in the position of the min- 
imum is supported by the Nielsen identity. As a function of the parameter 
£ and 0, the potential develops iso-potential curves, i.e. there is a function 
0(00)0 that fulfills 0(00; Co) — 0o an d according to (JHJ) 

e|v A e,- e |»V +t »V- ft »V +t Jv-0. ,46) 

An explicit expression for the function Co is given in ( I15p . The leading 
contribution obtained by evaluating the diagram shown in Fig. [1] reads 

1 2 t m f d3 P 



Co 2& J (2ir) 3 p 2 — m 2 p 2 — m FP ' 
that gives 

C^l. ?W ^ g ^T, (48) 

where the last expression is obtained when evaluating the masses close to the 
broken phase. This agrees nicely with our findings in ( H2|) and the relation 
C = ^d(j)/d^. Using ( HHj) . it is also possible to verify explicitly that the 
Nielsen identity for the effective potential ( 1T4|) is indeed satisfied also away 
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from the broken minimum, up to higher order corrections. Namely, taking 
the ^-derivative of the effective potential given in (I40p one finds 

*~dt = -^9 2 ^ 2 T(m x -m FP ). (49) 

The right-hand side is of higher order close to the minimum where m x m 
nipp, such that the value of the effective potential in the minimum is £- 
independent, up to higher order corrections, as discussed before. Using that 
1^- ~ 2(f)(m x — m 2 FP ) + 0(g 4 , Xg) in combination with (I48j) then shows that 
the Nielsen identity is respected also away from the broken minimum. 



6 Gauge Independence in the Gradient Ex- 
pansion 

Let us return to the relation (1251) that ensures the gauge independence of 
the effective action in the gradient expansion. An explicit calculation of the 
wave function correction appearing in the derivative expansion of the effective 
action yields (see Appendix [B]) 

ml 2 2 
1 1 

16m| m x + mp mpp + my 

where we have defined 

Within our counting, 8 is of order g and hence subleading. Still, we 
keep these corrections since unlike two-loop contributions they are linear 
in the temperature and (as we shall see) cancel among themselves. The £ 
dependence of Z arises via the masses m x and mpp, apart from the explicit 
dependence on £. Formally, the ^-dependent terms are of the same order as 
the ^-independent ones. In the limit £ — > 0, the wave function correction 
agrees with the one obtained in Landau gauge (see appendix A of [27]). 



9 2 T 11 



37r \ 8mp 



+ 1 



m x + mpp 
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Note that close to the broken minimum where m x m m FP , the leading 
^-dependent terms cancel and the £ dependence becomes suppressed. This 
behavior is precisely the one expected from the Nielsen identity fl25|) . because 
close to the broken minimum all terms on the right-hand side are at most of 
order g 2 . Namely, the term proportional to dCo/d(f> is of higher order because 
close to the broken minimum Cq becomes approximately <fr- independent. The 
term involving C dZ/d(j) is suppressed because both factors are of order g, 
and the term proportional to the derivative of the effective potential also 
has to vanish in the broken minimum by definition. Finally, the last term is 
suppressed as well because the Higgs mass is of order A ~ g 3 . In order to 
check explicitly that the Nielsen identity (I25p is also fulfilled for field con- 
figurations away from the broken minimum, we have computed the leading 
contributions to the coefficients D and D contributing to the derivative ex- 
pansion ff24l) of the Nielsen coefficient C(x) for the full effective action. It 
is possible to relate these coefficients to Feynman diagrams obtained from 
attaching one or two external Higgs field lines to the one-loop graph shown 
in Fig. [T] that corresponds to the leading contribution to C(x). The diagrams 
and the explicit expressions are shown in Appendix [Bl 

One can check that, with these coefficients (198]) . the result for Z from 
Eq. (}50l) . and Cq from Eq. (1481) . the Nielsen identity (|25|) for the correction of 
the kinetic term is indeed satisfied to order g 2 . For this calculation, we found 
it helpful to express also the leading contributions to the field derivatives of 
Cq and V in terms of the Goldstone and ghost masses, 

dCo „ g 2 Ttm 2 x -m 2 FP (l + S/0 
dcf> 8tt m x (m x + m FP ) 2 

and dV/d<f) ~ 2<f)(m x — m 2 FP ), as well as d 2 V/d<p 2 ~ 2(m^ — m FP ) + 4m FP 5/C, 
up to corrections of order 0(g A ,Xg). Two numerical examples are given in 

Fig. m 

In conclusion, we have the following findings: 

First: The function Cq is close to the broken phase of order g^f£T and 
hence the £ dependence of the critical field value <p c is subleading. Another 
important consequence is that the first-order nature of the phase transition 
is a gauge-independent statement [51 [6] in the regime where the perturbative 
expansion converges %. 

5 In our current scheme this also entails the condition that £ is not much larger than 
0(1). 
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g = 1/3 




Figure 7: The different contributions to the Nielsen identity (|25p for Z{4>) an d 
their sum. The labels (a) to (e) are the five contributions in (|25p from left to right. 
The parameters of the upper plot are g = 1/3 and A = 0.015, while the lower uses 
g = 1/10 and A = 4.05 x 10 -4 . This ensures 4> C /T C ~ 1 in both cases. The gauge 
parameter is £ = 0.5. 
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Second: The Nielsen identity (I25p for Z is in leading order (g 2 ) fulfilled 
everywhere. The derivative dC /d(j) is of order g 2 T/<p^. Hence, according 
to (j25p t^dZ/d^ can be at most of order g 2 if dZ/dcj) is at most of order g. 
From the explicit result (150]) it turns out that, in the broken phase, indeed 
C,dZ/dC, is of order g 2 while dZ/dcf) is of order g. Finally, D is of order 1/g, so 
all terms in f ]2"5"]) are of the same order and their sum to this order vanishes. 

Third: Even though we cannot check the equation (1251) in full general- 
ity to order g 3 , the one- loop terms are the only contributions to this order 
that are linear in the temperature. These contributions cancel among them- 
selves. The remaining terms in ( 125]) to order g 3 (denoted by E in Fig. [7]) 
result from the contributions proportional to T 2 coming from the product of 
one-loop terms in CodZ/dcf). These have to cancel against genuine two-loop 
contributions that we did not calculate. 

Fourth: Close to the symmetric phase the functions Z, D and D diverge 
what makes the gauge dependence introduced by subleading contributions 
large. However, this does not indicate a breakdown of perturbation theory. 
What breaks down is the gradient expansion of the effective action T and of 
the function C that was used to arrive at the relation (|2"5"|) . In principle, if the 
gradient expansion was avoided, explicitly gauge-independent results could 
be obtained. Nevertheless, in practice the gradient expansion is typically 
used to determine observables like the sphaleron energy or the tunneling 
rate perturbatively. In the next sections, we discuss to what extend this 
introduces a gauge dependence to the conventional analysis found in the 
literature. 



7 Thermal Tunneling 

We start the discussion of thermal tunneling by assuming that the gradi- 
ent expansion of the effective action is sound and only later investigate to 
what extent this really holds true. The rate of thermal tunneling is mostly 
determined by the action of the so-called bounce solution of the effective ac- 
tion [28, 29J. The (Euclidean) action with 0(3) symmetry is in leading order 
given by 

r = — J d PP 2 (<9>v>z(0) + v(<f>)) . (53) 

To simplify the discussion, we assume that the temperature is close to the 
critical temperature where the solutions of the equation of motion is in the 
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'thin-wall-regime'. In this case, the field only changes in a region p ~ R 3> L, 
where R denotes the size of the nucleated bubbles while L denotes the wall 
thickness of the bounce solution [29]. In this approximation, the action can 
be rewritten as 

Air f Air 
T = Y R 2 J dp {{d^f Z{<j>) + V{<j>j) - —R 3 e , (54) 

where e denotes the potential difference between the symmetric and the bro- 
ken phase. The first integral is invariant under a shift in p leading to a 
conservation of 

(d p (j)) 2 Z{<t>) - V{<t>) = const . (55) 

For the configuration that dominates the path integral in the tunneling pro- 
cess the constant on the right side vanishes. This gives for the action 

r = !VR 2 -— R 3 e, (56) 

where we defined the wall tension^ 

a = J dp ((d p <f>) 2 Z{(j>) + 1/(0)) = j ' dxj>s/V{(j>)Z{(j>) . (57) 

Finally, the bubble size R can be obtained by extremizing this expression, 
what gives the well-known result 



2(7 167TCT 3 . , 

R = —, r = -. 58 

This little exercise shows that a gauge-independent tunneling rate can in 
the thin-wall regime only be obtained if the wall tension is gauge-independent 
(the potential difference e is obviously gauge-independent due to Nielsen's 
identity). At first sight, it seems to be impossible that the wall tension is 
gauge- independent: In the effective potential the gauge dependence starts 
at relative order g, while in Z the gauge dependence first occurs at relative 
order g 2 . However this argument is not sound, since the gauge dependence 
arising from the potential is further suppressed in the wall tension. To see 



6 Our definition of the wall tension differs from the usual definition by a factor V2 due 
to our conventions for the kinetic term. 
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this, one can ignore the sub leading contributions in Z to order g 2 . This gives 
then for the wall tension 



A r<t>c J r&c J r<t>c An 

Z,i-a~ d4&y/V~ d^Coi-W-- d^^W. (59) 



d£ Jo d£ 

Here we used that the potential vanishes in both phases in the thin-wall 
approximation. dCo/dip is in the broken phase of order g 2 such that the 
gauge dependence of the effective potential can in the wall tension cancel 
against the gauge dependence in Z. 

Is this result specific to the thin-wall regime and how does this fit to- 
gether with the statement ff25|) ? On first sight, f )25|) only ensure the gauge 
independence of the action in the gradient expansion, or more specifically 
as long as Zd^(f)d^(j) <C V. On the other hand, for the tunneling bounce 
solution both terms are of the same order and this expansion does not apply. 
However, in the present context there is a way to derive a similar relation as 
( |25|) for the tunneling bounce without this constraint. Starting from ([8]), a 
functional derivative with respect to j(x) yields 

£^ = C(*,0(*U), (60) 

in case (p(x) fulfills the equation of motion, j(x) = 0. So, there is a class of 
solutions <£>(x,£) with 

^r(0(z),f)Ux)= V (z,e) = °- (6i) 



Using this relation in (I53p yields 
= AnT [ dpp 2 (dpcpdp 



^dC + dZ c + ^dZ 
dip dip d^ 



So qualitatively the same picture emerges as in the thin-wall regime. To 
relative order g, the shift ( 16 Up in combination with ( I14p ensures the gauge 
independence but only because dC/d<j) is of order g 2 while C is of order g. 

Using the expansion ([21]) in flBTj) then leads to a relation similar to (125 p . 
Unlike the condition Zd^d^tp V, the expansion (pMI) is well justified in 
the present context at least in the broken phase. Nominally the expansion 
parameter is p 2 /m 2 such that with p 2 4> 2 ~ V ~ g 3 4>i and masses of order 
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m 2 ~ 2£g 2 (J) 2 the gradient expansion of ( I24p is valid for 3> \fg]~£,4>c- I n 
particular, the gradient expansion gets worse for small £. 

Now notice that in the last section we have shown the relation £dV/d$, = 
CodV/dcj) only up to terms of order g 5 and that dD/dcf) is like dCo/d<p of 
order g 2 . Hence to order g 5 , the relation (ISTj) is not exactly equivalent to 
f l2"5"j) . In any case, the range of validity is not the same for both relations, 
since ( l2~5"j) is valid for small gradients, while in the derivation of (ISTj) we used 
the equation of motion. 

From above discussion, it seems that the gauge dependence of the wall 
tension is suppressed at least by g 2 if Z is ignored and that the inclusion of 
the Z factor might even postpone it to order g 3 . However, the breakdown of 
the gradient expansion in the symmetric phase prevents a gauge-independent 
determination of the wall tension in (I5TJ) to this order. In practice, the ex- 
pression we obtained for the wave function normalization Z can even become 
negative for <fi ~ g<p c such that no reasonable results can be obtained if the 
leading non-trivial order of Z is taken into account. Also dC^/dcf) is of order 
y/g close to the symmetric phase what indicates that a large gauge depen- 
dence arises if the gradient expansion is employed. 

The gauge dependence arises mostly from the symmetric phase and inte- 
gration of Eq. ( 159]) leads to the estimate 

for £ not too small. Here m x denotes the Goldstone mass for = 0, that 
coincides with the Higgs mass rrih in the symmetric phase. Compared with 
the leading order result 

a ~ -m x 4> 2 c , (64) 
the uncertainty from the gauge dependence scales as yf\ ~ g 3 ^ 2 , 

(65) 



£ da 1 m y T c 



ad^ 8tt <p 2 c 
This leads to a logarithmic dependence on £, 

(7(0/(7(1) - W/ 2 ln(£), (66) 
which nicely fits with our numerical findings shown in Fig. [HJ Especially, 
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Figure 8: Dependence of the wall tension a on the gauge parameter £. The 
three values ao, oz and an denote the wall tension deduced, from V only, from V 
and the naive Z and from the full momentum dependence in II, respectively. For 
comparison, also the critical Temperature T c and the critical field expectation value 
4> c are shown. The residual gauge dependence of the wall tension is of relative order 
g 3 / 2 ln(£), while for T c it is very small, of order g 3 . In contrast, the £ dependence of 
4> c is expected from the Nielsen identity, see Eq. ([4~2jh and its leading contribution 
scales as gy/£. The parameters of the upper plot are g = 1/3 and A = 0.015, while 
the lower uses g = 1/10 and A = 4.05 x 10~ 4 . 
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this means that for very small values of £ the gauge dependence becomes en- 
hanced. For £ < g, the logarithm is ultimately cut off, and the £ dependence 
is of the relative order g 3 ! 2 ln(g). Therefore, it seems that a choice £ 3> 0(g) 
is slightly preferable for the computation of the wall tension compared to e.g. 
Landau gauge. 

Including the wave function correction Z to the kinetic term in the deter- 
mination of the wall tension mainly leads to a gauge-independent shift. This 
shift can be estimated as (details are again given in Appendix [C]) 

1 1^/2 

This scales as Aa/cx ~ g what is parametrically larger than the gauge- 
dependent contributions to the wall tension that scale as g 3 ^ 2 . The gauge 
dependence of the wall tension is not improved systematically by the inclu- 
sion of Z. Especially, it is still of order g 3 ^ 2 . Nevertheless it turns out that 
a cancellation of the ^-dependent contributions occurs for £ of order one. 
These findings are supported by our numerical results that are presented in 
Table [1] and Fig. [HJ Some analytical details of this discussion are given in 
Appendix O 

At this point, we would like to comment on [31] that compared the wall 
tension in the Landau and Feynman gauges for an effective three dimensional 
theory. There it was also found that including the Z factor does not lead to 
a systematic reduction of the gauge dependence. Nevertheless, the analysis 
showed that the gauge-independent contributions from Z where as large as 
the gauge- dependent ones. This discrepancy compared to our analysis is due 
to the fact that the Goldstone mass was not resummed what leads to a larger 
gauge dependence in the effective potential as stressed in section 

The fact that the inclusion of Z does not postpone the gauge dependence 
of the wall tension to the relative order g 2 or even g 3 is related to the break- 
down of the gradient expansion in the symmetric phase, which leads to a 
divergence of Z for — > 0. To avoid this problem, we would like to check the 
gauge independence of the effective action in the symmetric phase without 
resorting to the gradient expansion. In order to achieve that, one can expand 
the effective action in <f> around the symmetric phase and obtains 

V = T S ' (68) 
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with 



n(p) =p 2 + m 2 h + n 2 (p). 



(69) 



The gauge independence of the effective action (fl2|) then implies in leading 
order 

e-^-^^-ga, (70) 

where C (p) is now understood to be expanded in instead of p. In the limit 
— >■ 0, the only contribution at one-loop order to II 2 is the last diagram 
depicted in Fig. [12] involving the gauge and Goldstone boson and several 
tadpole diagrams that however have no momentum dependence. Explicit 
calculation shows 



cffl 2 



g 2 n 

4tt 



-m y + 



2(p 2 + m 2 
P 



-arctan^/ra^ 



(71) 



and dC / dip in accordance with (1701) . Notice also that II 2 is finite in the limit 
p — » if the limit — > is taken first 



dlh 

d£ 



g 2 n 

47T 



m x + 



8 p z 
3 m 



x J 



(72) 



The mass term agrees hereby with the one derived from the potential 
while the kinetic term does not 



^symm 1 



dp 2 



1 



p=0 



9 2 T 



(2 + 



(73) 



Since the potential agrees independent from what limit (0 — > or p — > 0) is 
taken first, it is tempting to reduce the gauge dependence of the wall tension 
by interpolating between the two different kinetic terms Z and Z sjmm . In 
order to implement this idea, we replace the wave function correction in 
Eq. (EZD by 

Z(M = (74) 

Here n 2 (0,p) is the full one- loop self-energy for a general expectation value 
of the Higgs field, which approaches the self-energy II 2 (p) discussed above 
for — y 0. On the other hand, Z(<ft,p — )■ 0) agrees with the correction to 
the kinetic term from Eq. ( l50|) . For any non-zero value of the momentum p, 
Z(<f),p) has a regular behavior for — )• 0, and indeed interpolates between Z 
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(7 = 1/3 


(7 = 2/3 


A 


0.0075 


0.015 


0.025 


0.06 


0.1 


0.15 


(£=1) 


6.29 


1.78 


0.74 


149 


31.2 


12.9 


5 (T 


0.11% 


0.33% 


0.79% 


0.33% 


0.70% 


1.4% 


°z (£ = 1) 


6.16 


1.70 


0.68 


143 


28.8 


11.3 


8a z 


-0.10% 


-0.25% 


-0.56% 


-0.21% 


-0.46% 


-0.97% 


M£ = i) 


6.16 


1.70 


0.69 


143 


28.8 


11.3 


5<Jn 


-0.07% 


-0.16% 


-0.32% 


-0.13% 


-0.25% 


-0.52% 



Table 1: Numerical results for the wall tension in units of 10 _3 t> 3 , and the 
shift 8 a for a change of the gauge parameter £ from 1 to 0.1. The rows a , 
<y z and an denote the wall tension deduced, from V only, from V and the 
naive Z and from the full momentum dependence in II, respectively. 

and Z symm as a function of <fi provided that the momentum is chosen small 
enough, p < m^. In Tabled] we show the resulting expression for the wall ten- 
sion, where we used p = m^/fiir) as the momentum cut-off. Unfortunately, 
we found that qualitatively the gauge dependence did not improve signifi- 
cantly by doing so, although it is slightly reduced on a quantitative level. 
A complete cancellation of the gauge dependence at the relative g 3 ^ 2 level 
seems to require the use of the full momentum dependence of the effective 
action without resorting to the gradient expansion. 

8 Sphaleron numerics 

In this section we briefly discuss the sphaleron energy following [52]. The 
sphaleron is a static Higgs-gauge configuration that is a saddle point of the 
action (which reflects the energy of the configuration). It has Chern-Simons 
number | and is situated half-way between two gauge vacua. Our toy model 
does not contain a 577(2) gauge sector and hence no sphaleron transitions, 
but the Higgs potential in the Standard Model has essentially the same fea- 
tures of the Abelian Higgs and it is reasonable to feed the potential (HOI into 
the equations of motion of the sphaleron to estimate the gauge dependence of 
the sphaleron in the Standard Model. In the conventional analysis, the main 
difference between the Abelian and the non-Abelian model in terms of the 
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strength of the phase transition is that there are three times as many gauge 
bosons (and ghosts and Goldstones) contributing to the cubic term hence 
strengthening the phase transition. We will mimic that by also presenting 
results for larger than observed gauge couplings. 

In the non-Abelian case, perturbation theory is plagued by Linde's prob- 
lem in the symmetric phase such that the expansion in the coupling constant 
becomes questionable for small Higgs vevs. Besides, compared to the tun- 
neling rate discussed in the last section, the convergence of the gradient 
expansion is even more problematic. The gradient expansion is formally an 
expansion in the parameter p 2 /m 2 where the relevant mass is the one of the 
gauge bosons. While for the tunneling bounce this is 0(g), in the case of 
the sphaleron the expansion parameter is 0(1) and hence not suppressed by 
any coupling constant even in the broken phase. However, numerically the 
coefficients D and Z — 1 are slightly smaller than Cq and V e ff (in units of 
m 2 A 4> 2 c ) such that higher orders can be neglected if this trend continues. This 
issue is to certain extend unrelated to the gauge dependence. 

At the same time, the sphaleron energy is proportional to <p c such that 
the gauge dependence stemming from <p c cannot possibly be canceled by 
the gauge dependence of the effective action solely in the surrounding of 
the broken phase. In order to quantify our lack of knowledge on Z and 
V e ff in the symmetric phase, our strategy is to set Z to 1 in the numerical 
analysis and use the gauge dependence of V e ff to estimate the impact of 
those contributions to the sphaleron energy. 

The differential equations to solve when a spherical Ansatz is used read 

C 2 = 2/(1 -/)(!- 2/) -I C 2 ^ 2 (1-/), (75) 
d 2 dh 2 1 dV h 

dC C Tc = 2h{1 - f) + ?^' (76) 

with the asymptotic behavior 

f^a( 2 , h^P( for C^O, 

and 

/->l- 7 exp(-C/2), h^l- 6 -exp(-K() for C ^ oo , (77) 

where the parameter k is given by k 2 = and the rescaled potential is 
defined as Vh(h) = V e ff(h ■ <p c ) while the rescaled coordinate is ( = g<p c \x\. 
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= 1/3 


= 2/3 


A 


0.0075 


0.015 


0.025 


0.06 


0.1 


0.15 




0.965 


0.549 


0.325 


2.07 


1.03 


0.535 


Tjv 


0.454 


0.575 


0.721 


1.00 


0.94 


1.05 


8<p c 


0.3% 


0.7% 


1.5% 


0.6% 


1.3% 


3.0% 




0.01% 


0.01% 


0.02% 


0.05% 


0.07% 


0.10% 



Table 2: Numerical results for the sphaleron energy. The row 5<p c contains 
the shift for a change of the gauge parameter £ from 1 to 0.1. The row 
contains the corresponding change in sphaleron energy due to a change in 
the shape of the Higgs potential. 



We solve the equations numerically with a shooting algorithm similar 
to what is used to find the bounce solution of the tunneling action. The 
parameters a and /3 are chosen and the equations are solved from some 
position ( e close to the origin to a value £ w ~ 0(10). The parameters a and 
(3 are then varied and we search for simultaneous zeros in the functions 



D f = 1-/ + 2/', 

D h = (K(l-h) + ((h)'. 



(78) 



These two conditions ensure that the numerical solutions smoothly match 
the asymptotic behavior given in (ITT)) . The sphaleron energy is then given 
by 



E 



4\Z2~7T0 C 

9 



4(f) +£(/(!-/)) 



+ 1 2 (2 (tc I +(MW)) 2 + CW, 



(79) 



Some numerical results are given in Table |2j Main impact of the sphaleron 
energy has the gauge-dependent shift in <f) c . The gauge dependence in the 
critical temperature T c and the shape of the effective potential are sublead- 
ing. Moreover, close to borderline case of sphaleron washout C /T C ~ 1, the 
uncertainty in the sphaleron energy never exceeds a few percent for Standard 
Model values of the gauge coupling. In summary, the uncertainty stemming 
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from the residual gauge dependence is subleading compared to corrections 
coming from two- loop contributions to the effective potential [331 HO] • This 
gauge dependence is inherited from the critical vev 



and scales as g. To remove this gauge dependence one would probably need 
to include the next-to-leading order of the kinetic term of the gauge-bosons. 
There, a sizable gauge dependence is expected in order to ensure a gauge- 
independent position of the pole in the gauge-boson propagator. Besides, the 
breakdown of the gradient expansion should lead to even more severe effects 
than in the case of thermal tunneling as discussed above. This will further 
complicate the determination of the sphaleron energy with an accuracy be- 
yond the bound (IHUj) . 

9 Summary 

Let us summarize our findings concerning the gauge dependence of the effec- 
tive action in the Abelian Higgs model in i?^-gauges. We explicitly demon- 
strated various Nielsen identities in the regime where the use of perturbation 
theory and the gradient expansion of the effective action is feasible. 

In particular, we have shown that the position of the minimum of the ef- 
fective potential transforms in leading order according to f|T4l) under a change 
in the gauge fixing parameter £. We would like to emphasize that this result 
could only be obtained by calculating the effective potential consistently to 
order g 3 using the counting g 3 ~ A. In particular, this required the resum- 
mation of contributions to the Goldstone boson masa3 of order m^T (where 
TiiA denotes collectively the different masses of the gauge bosons). 

Furthermore, we have demonstrated that the off-shell effective action in 
gradient expansion transforms according to the Nielsen identity (I25I) . How- 
ever, this relation cannot guarantee the gauge independence of vacuum tran- 
sitions for several reasons. First, the gradient expansion is not well justified 
in these cases since the contribution of the kinetic term to the action is of 
equal size as (or even larger than) the contribution from the scalar potential. 

7 At this point of the analysis, we depart from [MJ 133 IM] that argue for a large gauge 
dependence of the effective action. 
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Second, even for small gradients the gradient expansion and the relation (I25p 
break down at some point in the symmetric phase and vacuum transitions are 
also sensitive to this regime. Compared to the analysis at zero temperature 
[TT] . notice also the additional contribution involving D that was missing but 
also not important in the analysis presented there. 

Finally, we discussed the gauge dependence of the tunneling action. Using 
the established procedure to calculate the tunneling action perturbatively 
(meaning a canonical kinetic term and the appropriate effective potential), 
we found for the gauge dependence of the wall tension the estimate 

ad£ 8tt 02 ' \ ) 

what scales as \f\ ~ g 3 ^ 2 (the Goldstone mass m x is evaluated in the symmet- 
ric phase). Including corrections to the kinetic term Z leads to corrections of 
the same order but the gauge dependence is not persistently reduced. This 
is due to the fact that the latter corrections are sensitive to the effective ac- 
tion very close to the symmetric phase where the gradient expansion breaks 
down. 

Since we did not arrive at an explicitly gauge-independent result for a, 
this leaves the question what is the best gauge to chose. Our results develop 
the strongest £ dependence for £ ~ 0. For £ > g the wall tension computed 
using a canonical kinetic term depends logarithmically on £, while when 
including the correction to the kinetic term the £ dependence partially cancels 
and a is rather insensitive to £ for £ ~ 1. This could indicate that a £ value 
of this order is the appropriate choice and not Landau gauge (that is mostly 
used in the literature and is reproduced in our case by £ — > 0). Nevertheless, 
quantitatively the dependence on £ is rather small and including the two-loop 
contributions to the effective potential will probably have a larger impact in 
most models (e.g. contributions from the gluons [3?]). We also identified 
gauge-independent corrections to the wall tension that arise from the kinetic 
term Z and scale as g. Also these contributions are more important than the 
gauge-dependent ones. 

Several non-perturbative studies of tunneling in the Standard Model were 
presented in refs. [201 EH] an d a detailed comparison with the perturbative re- 
sults in Landau gauge was given in [39] . Also there it was concluded that the 
corrections to the kinetic term and two-loop contributions to the effective po- 
tential are important to achieve a good agreement with the non-perturbative 



34 



results for the wall tension. Numerically, these corrections are far more im- 
portant than the gauge fixing dependence we discussed here. 

In case of the sphaleron, the convergence of the gradient expansion is also 
problematic. As a naive estimate of the gauge dependence of the sphaleron 
energy, we obtained 

I d_E ^ j^_Tc 

E di 167TV/2 0C' 

what is typically of a few percent. Over all, the effective potential enters in 
the sphaleron energy mostly via the position of its minimum. Hence, im- 
proving the gauge dependence of the sphaleron energy will probably require 
not only to go beyond the gradient expansion but also to calculate the gauge 
dependence of the wave function corrections of the gauge fields, what we did 
not attempt here. 

The estimates ( 1821) and ( 1671) are solely based on the Nielsen identity for 
the Higgs vev, C = ^dcf)/d^. Therefore, the results can be readily carried 
over to the Standard Model and some of its extensions as the two-Higgs 
doublet model or singlet extensions. Compared to the Abelian Higgs model, 
the electroweak sector of these models gives in leading order a gauge fixing 
dependence of the Higgs vev that is larger by a factor three. 

In conclusion, determining vacuum transitions in an explicitly gauge- 
independent fashion is mostly hindered by the breakdown of the gradient 
expansion of the effective action (in particular in the vicinity of the symmet- 
ric phase). Still, in the cosmologically most interesting regime with g <C 1 and 
4> C /T C > 1, the gauge dependence of the tunneling action and the sphaleron 
energy is rather small. The situation further improves when the essential cu- 
bic contributions to the effective potential do not solely arise from the gauge 
bosons (as in the present Abelian toy model) but e.g. from additional degrees 
of freedom (as in the light stop scenario [10| W\\ W2\ ) or from the tree level 
dynamics of an extended scalar sector |43j. 
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Figure 9: Feynman rules for scalars in the Abelian Higgs model. 

A Feynman Rules 

The Feynman rules for the scalar particles are shown in Fig. [9] while the 
Feynman rules involving the vector particles are given in Fig. [TUJ 



B Wave function corrections 

In this section, we present some results on the coefficient of the kinetic term 
Z in the i?^-gauge and also for the functions D and D. As mentioned in 
section HI derivatives with respect to the background field <fi are in the R^- 
gauge not related to diagrams with external Higgs fields. In order to remedy 
this issue, we replace in the gauge fixing terms the background field </> by 
4>+h/ \/2 where h will be treated as an external field. The function Z is then 
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Figure 10: Feynman rules for vectors in the Abelian Higgs model. 

related to the two-point functions involving h and h via 

d 1 



7/1U 



dp^dp v 



(n hh + u hh + u- hh + r%) . (83) 



The additional Feynman rules involving the external field h are shown in 
Fig. [TTJ For the sphaleron and thermal tunneling, we only need the spatial 
components Z n at finite temperature. 

On the one loop level the possible diagrams are of the form depicted in 
Fig. [12] where the external lines can be either h or h. On one loop level it 
is more practical to first perform the sum in f )83|) before the integrals are 
evaluated. The momentum dependence stems from the integrands 

1 1 



(I + p/2) 2 + m 2 (l- p/2) 2 + m 2 

1 1 p 2 (p- I f 



(l 2 + m 2 ) 2 2(/ 2 + m 2 ) 3 (/ 2 + m 2 ) 4 
1 1 p 2 1 p 2 m 2 



(l 2 + m 2 ) 2 6(/ 2 + m 2 ) 3 3(/ 2 + m : 



2U 



(84) 
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Figure 11: Additional Feynman rules for the kinetic term in the effective action. 
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Figure 12: The diagrams contributing to the kinetic term of the effective action, 
that can be evaluated using 

L B (y) = 7^J2j dx 
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and a similar expression for the the integrals with higher powers. Notice that 
only the zero mode n = contributes to the leading term and we neglect 
higher Matsubara modes in the following whenever this holds true. 

We first restrict ourselves to the broken phase where the masses are given 

as 

m 2 A = 2g 2 <p 2 = m 2 = yT 2 , m 2 = m 2 FP = £m 2 . (86) 
Then one finds for the ghost and Goldstone loops the contributions 

^ff>ey(L B (£y)+2ZyM B {Sy)). (87) 

The diagram involving two gauge fields requires the evaluation of products 
of the polarization projections P Miy (/) = IH V jl 2 

P^(l+p/2)P, u (l-p/2) ~ i - ^ + + (5(p») 

- l ~ll? Jr0{p " ) - (88) 

and again the expansion of the integrand (151)1 . The novel integrals are given 
in appendix [D) The contribution from the term involving the p M polarization 
twice gives 

±g 2 ey (-L B (£y) - 2^yM B ^y) -4L B {Zy, £y, 0)) . (89) 
The contribution involving one p M polarization yields 

~g 2 tyL B (ty,y,0), (90) 
and the contribution with no p M polarization is 

2g 2 y (-~L B (y) - yM B (y) - ~L B (y,y,0)\ . (91) 
Finally, the contributions from the mixed gauge-Goldstone loop give 

4/eQ^(&/,&/)J , (92) 

and 

^9 2 (^K B ^y,y)^ • (93) 
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In the broken phase, where m 2 FP = m 2 the contributions proportional to 
Lb in dHZ]) and (I89p cancel each other. The third term involving Lb in (l89p 
cancels in leading order against f )92|) . Besides, the terms in (190]) and f )93|) 
combine to 

^ 2 (^L B (^,y,0) + ^(^,2/)) -3^. ( 94 ) 

what is gauge-independent. Together with the contribution in ( 19T|) this yields 
the final result 

Z ~ 1 ^— . (95) 

Taking in addition the resummed gauge boson propagator 



- ml ^ p 2 -m 2 T ^ % 2 - m 2 FP ^ ' 1 J 

into account, where P°„ = p^p u /p 2 , P^ u = - u^u u - (p^p^/p 2 , P^ v = 
9»u - (p^Pu)/p 2 ~ PjL, P^ = P^~ u^up), = (1, 0, 0, 0), m\ = m\ + ag 2 T 2 } 
and m\ = 2g 2 <f) 2 , we find 

9 2 (97) 



The above calculation can be extended to a configuration in field space 
away from the broken minimum. In that case one cannot set the Goldstone 
mass m x and the mass m FP of the ghost and of the time-like gauge boson 
polarization to be equal. In that case we obtain the result given in Eq. ( I5U]) 

A similar calculation can be done for the factors D and D in (124"]) . These 
are obtained from calculating the diagrams shown in Fig. [13] and Fig. HU re- 
spectively, expanding in powers of the external momentum p M , and extracting 
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Figure 13: The diagrams contributing to the function D. 



the contribution proportional to p 2 . This yields the result 
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Figure 14: The diagrams contributing to the function D. 

C Wall tension 

In this section we discuss the gauge dependence of the wall tension in more 
detail. 

First, we show that an explicit calculation reproduces the dependence 
advocated in f )63|) by use of the Nielsen identity. The wall tension (without 
including the Z factor) is given by 

a = J dtj>y/Vv + 5V , (99) 

where we split the potential into a gauge-independent piece V and the con- 
tribution arising from the Goldstone bosons ( I4ip . The contribution to the 
wall tension from the broken phase is nearly gauge-independent as long as 
the integration boundaries are adapted consistently. Close to the symmetric 
phase, two sources of gauge dependence arises. The first stems from the up- 
per boundary of integration and the second from SV. Due to 5V Vo one 
finds 



1 Jo VV 



(100) 
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and (j) is arbitrary as long as 2g 2 £0 2 3> m 2 . Let us first evaluate this integral 
in the regime <C <p c . This requires £ 3> A/2g 2 and leads to ~ m 2 2 and 
m 2 = A0 2 /4 = const. The corresponding integral can be evaluated using 

f s Ht 3 4 

/ ^±((l + x 2 ) 3/2_ x 3_ 1 ) = =_ log (^_ +log2 + 0(bg(5)/x). (101) 

Jq x 2 6 

We are mostly interested in the two leading terms. The constant term 
do not lead to a £ dependence and the subleading corrections are of order 
A log(g)/g 2 £ ~ glog(g)/£. In the wall tension, one cannot chose paramet- 
rically smaller than <f) c , but one can neglect the range of integration from 
zero to a few, since it does not contribute to the two leading terms. This 
makes the result meaningful, since the effective action cannot be trusted for 
x < 1 due to the breakdown of the gradient expansion. This also allows one 
to expand the numerator of the integrand and to include the full mean field 
potential and not just its linearization in the symmetric phase. The arising 
integral is 

I , dX / A -x(l - x/x e )(l - 2x/x c ) - 1 ) 
J x(l-x/x e ) \2 1 ' cA 1 c) ) 

3 

= -x(l - x/x c ) - \og(x) +log(l - x/x c ) 

+ const + C(log(x) /x) . (102) 

After this procedure, x can be chosen as a fixed multiple of x c , e.g. x = 
x c /2 = \j2g 2 £ i /\. The linear term cancels then against the contribution 
from the boundary of integration (what can be checked using f|48|) ) while the 
logarithmic term reproduces the relation (I63|) 

Aa ~ T^T c m 2 log(£) . (103) 

487T A 

The constant term is gauge-independent and corrections to this relation are 
of relative order glog(g)/£. 

Next, we examine if the logarithmic terms cancels against the contribu- 
tions arising from the wave function corrections Z. For larger values of x, 
the integrand can be expand as 

Aa z = - [ d(f)SZVV. (104) 
2 Jo 
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This integral still diverges, but one can extract a gauge-independent contri- 
bution 



3n \8rriT 16m L 



>^l^-51' (105) 



We will deal with this contribution later. 

Since the leading contributions from the remaining SZ is of order log(x), 
the potential can be linearized from the start and the integration of the 
^-dependent terms of 5Z in (|50|) leads for large x to a contribution 



* " 12^™* {--2 + —Z + (T^ ) ■ ( ° 6) 

Note that the right-hand side is regular for £ > 0, especially it has no pole 
for £ = 1. It is also possible to estimate the £ dependence by evaluating 
the ^-derivative similarly to Eq. (|59|) . but including Z and using the Nielsen 
identity Eq. ( )25|) . This leads to the same result as shown above. For £ not 
too small, the £ dependence of Ado and Aaz are of the same order and have 
opposite signs. However, there is no systematic cancellation between the £ 
dependences of these two quantities. In any case, it turns out that it is not 
justified to neglect the integration in x from zero to a few, since the corre- 
sponding contribution is potentially as large as the one we just presented. 
This contribution cannot reliably determined due to the breakdown of the 
gradient expansion close to the symmetric phase. Ultimately this prevents 
us from obtaining a gauge-independent result for the wall tension. 

Finally, consider the gauge- independent piece f ll05p . The corresponding 
contribution to the wall tension can be evaluated using 




ax<2 I n ^ / i«for/3 = /, n7 s 

" ( x 2 + ^3/2 ~ J x ^ ~ x >- |o.0146 a for ~ 1 ' ^ } 

Hence, the longitudinal gauge bosons contribute via the wave function cor- 
rection a term 

Aa~-^^gm x T c c , (108) 

to the wall tension. Notice that this is of order g 5 ^ 2 and parametrically larger 
than the gauge- dependent contributions we estimated before. 
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D Integrals 

In the following we list the used one-loop integrals. 



My) = Yl J dx x2 log ( 4ff2n2 + x2 + v) 



const + -^y - ^y 3/2 + 0{y 2 logy) , (109) 



My) = ^E/^y 



47r 2 n 2 + x 2 + y 

A/ B ( y )~-i-yV2 + Oology), (110) 



1 ' +0(\ogy), (111) 



I f 1 

L B (y 1 ,y 2 ,y 3 ) = ^ / dx x 2 



x 2 + y l x 2 + y 2 x 2 + y 3 



1 



L B (y) = L B (y,y,y), 
M Biv) = ^y^^^.^ + o^). (U2) 
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